CGNS API和FVM在非结构混合网格计算中的应用

用CGNS API（CFD General Notation System Application Programming Interface）作为非结构混合网格求解器的前处理和后处理，用FVM（Finite Volume Method）作为偏微分方程求解方法．在前处理过程中，用hash表法对内部网格面和边界网格面进行编号，并计算出相应的几何信息，以满足FVM求解器的需要．从FVM求解器计算出来的各种场信息可以写入原来的CGNS文件，该文件可以被许多专业商业后处理软件（如Tecplot，Fluent，CFX等）读取和进行可视化；对于求解器，用基于网格中心的FVM及SIMPLEC（Semi Implicit Methodfor Pressure Linked Equation Consistent）方法求解压力速度耦合．最后给出两个说明算例．     

非结构网格上温度扩散方程的能流计算方法

讨论非结构网格上温度扩散方程的能流计算方法.应用有限点方法(Finite Point Method,简称FPM)导出基于有限点两点公式和三点公式的能流计算公式,该公式适用于任意多边形及非匹配网格等非结构网格;给出网格角点温度新的计算公式.数值试验表明:基于两点公式的离散解和基于三点公式的离散解均具有平方阶的收敛速度;基于三点公式的离散解的精度总优于基于两点公式的离散解.     

A conservative unstructured scheme for rapidly varied flows

This article introduces a new semi‐implicit, staggered finite volume scheme on unstructured meshes for the modelling of rapidly varied shallow water flows. Rapidly varied flows occur in the inundation of dry land during flooding situations. They typically involve bores and hydraulic jumps after obstacles such as road banks. Near such sudden flow transitions, the grid resolution is often low compared with the gradients of the bathymetry. Locally the hydrostatic pressure assumption may become invalid. In these situations, it is crucial to apply the correct conservation properties to obtain accurate results. An important feature of this scheme is therefore its ability to conserve momentum locally or, by choice, preserve constant energy head along a streamline. This is achieved using a special interpolation method and control volumes for momentum. The efficiency of inundation calculations with locally very high velocities, and in the case of unstructured meshes locally very small grid distances, is severely hampered by the Courant condition. This article provides a solution in the form of a locally implicit time integration for the advective terms that allows for an explicit calculation in most of the domain, while maintaining unconditional stability by implicit calculations only where necessary. The complex geometry of flooded urban areas asks for the flexibility of unstructured meshes. The efficient calculation of the pressure gradient in this, and other semi‐implicit staggered schemes, requires, however, an orthogonality condition to be put on the grid. In this article a simple method is introduced to generate unstructured hybrid meshes that fulfil this requirement. Copyright © 2008 John Wiley & Sons, Ltd.     

Two splitting schemes for nonstationary convection-diffusion problems on tetrahedral meshes

Two splitting schemes are proposed for the numerical solution of three-dimensional nonstationary convection-diffusion problems on unstructured meshes in the case of a full diffusion tensor. An advantage of the first scheme is that splitting is generated by the properties of the approximation spaces and does not reduce the order of accuracy. An advantage of the second scheme is that the resulting numerical solutions are nonnegative. A numerical study is conducted to compare the splitting schemes with classical methods, such as finite elements and mixed finite elements. The numerical results show that the splitting schemes are characterized by low dissipation, high-order accuracy, and versatility.     

Evaluation of adaptive mesh refinement and coarsening for the computation of compressible flows on unstructured meshes

The compressible gas flows of interest to aerospace applications often involve situations where shock and expansion waves are present. Decreasing the characteristic dimension of the computational cells in the vicinity of shock waves improves the quality of the computed flows. This reduction in size may be accomplished by the use of mesh adaption procedures. In this paper an analysis is presented of an adaptive mesh scheme developed for an unstructured mesh finite volume upwind computer code. This scheme is tailored to refine or coarsen the computational mesh where gradients of the flow properties are respectively high or low. The refinement and coarsening procedures are applied to the classical gas dynamic problems of the stabilization of shock waves by solid bodies. In particular, situations where oblique shock waves interact with an expansion fan and where bow shocks arise around solid bodies are considered. The effectiveness of the scheme in reducing the computational time, while increasing the solution accuracy, is assessed. It is shown that the refinement procedure alone leads to a number of computational cells which is 20% larger than when alternate passes of refinement and coarsening are used. Accordingly, a reduction of computational time of the same order of magnitude is obtained. Copyright © 2005 John Wiley & Sons, Ltd.     

美式期权定价问题的变网格差分方法

本文提出一种求解美式期权定价自由边值问题的变网格差分方法.通过建立一个自由边界所满足的方程,利用变网格技术可同时求出期权的差分解和最佳执行边界.本文分别讨论了显式和隐式变网格差分格式,并给出了差分解的收敛性和稳定性分析.数值实验表明本文算法是一个非常有效的期权定价算法.     

High-Order Leap-Frog Based Discontinuous Galerkin Method for the Time-Domain Maxwell Equations on Non-Conforming Simplicial Meshes

A high-order leap-flog based non-dissipative discontinuous Galerkin timedomain method for solving Maxwell's equations is introduced and analyzed. The proposed method combines a centered approximation for the evaluation of fluxes at the interface between neighboring elements, with a Nth-order leap-frog time scheme. Moreover, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable under some CFL-like condition on the time step. The convergence of the semi-discrete approximation to Maxwell's equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with highorder elements show the potential of the method.     

A compact high‐order HLL solver for nonlinear hyperbolic systems

We present a new finite‐volume method for calculating complex flows on non‐uniform meshes. This method is designed to be highly compact and to accurately capture all discontinuities that may arise within the solution of a nonlinear hyperbolic system. In the first step, we devise a fourth‐degree Hermite polynomial to interpolate the solution. The coefficients defining this polynomial are calculated by using a least‐square method. To introduce monotonicity conditions within the procedure, two constraints are added into the least‐square system. Those constraints are derived by locally matching the high‐order Hermite polynomial with a low‐order TVD polynomial. To emulate these constraints only in regions of discontinuities, data‐depending weights are defined; these weights are based upon normalized indicators of smoothness of the solution and are parameterized by an O(1) quantity. The reconstruction so generated is highly compact and is fifth‐order accurate when the solution is smooth; this reconstruction becomes first order in regions of discontinuities. In the second step, this reconstruction is inserted in an HLL approximate Riemann solver. This solver is designed to correctly capture all discontinuities that may arise into the solution. To this aim, we introduce the contribution of a possible contact discontinuity into the HLL Riemann solver. Thus, a spatially fifth‐order non‐oscillatory method is generated. This method evolves in time the solution and its first derivative. In a one‐dimensional context, a linear spectral analysis and extensive numerical experiments make it possible to assess the robustness and the advantages of the method in computing multi‐scale problems with embedded discontinuities. Copyright © 2008 John Wiley & Sons, Ltd.     

Stokes问题各向异性网格Q2-P1混合元超收敛分析

讨论Stokes问题在各向异性冈格下的Q2-P1混合有限元方法，利用积分恒等式技巧得到了与传统方法相同的超逼近性质，同时基于插值后处理的技巧，构造了速度和压力的一对插值后处理算子，并且前者具有备向异性特征，从而导出了整体超收敛结果．